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Estuarine  Sediment  Budgets  for 
Chesapeake  Bay  Tributaries 

by  Julie  D.  Herman 


PURPOSE.  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  summarizes 
estuarine  suspended  sediment  budgets  developed  for  several  estuaries  of  the  Chesapeake  Bay  in 
Maryland  and  Virginia  (Herman  and  Eriedrichs  2010).  These  sediment  budgets  provide  better 
understanding  of  sediment  transport  patterns  and  pathways,  and  will  aid  in  improving  Regional 
Sediment  Management  (RSM)  practices  in  the  US  Army  Engineer  Districts,  Baltimore  and 
Norfolk.  Sediment  budgets  were  calculated  for  the  Patuxent  River,  Maryland;  York  River, 
Virginia;  and  a  partial  budget  was  developed  for  the  Potomac  River,  Maryland. 

BACKGROUND. 

Sediment  Budgets.  Sediment  budgets  account  for  the  sources  (areas  of  erosion)  and  sinks 
(areas  of  deposition)  of  sediment  in  a  specified  area  over  a  specified  time  frame.  The  components 
of  a  sediment  budget  may  vary  depending  on  spatial  and  temporal  scales,  local  variables,  and 
available  data.  Sediment  budgets  can  identify  types  and  locations  of  missing  data  that  could  help 
direct  future  research  needs  for  additional  data  collection.  Also,  sediment  budgets  may  assist 
significantly  in  improving  management  decisions. 

The  sources  and  sinks  in  a  sediment  budget  are  calculated  as  sediment  loads,  and  the  loads  are 
related  in  a  mass  balance  equation.  Eoad  is  mass  per  unit  of  time.  This  study  used  metric  tons  per 
year  (Mt/yr),  where  a  metric  ton  is  1,000  kg.  Eigure  1  displays  the  conceptual  model  of  the 
sediment  budget  for  Chesapeake  Bay  estuaries. 

Study  Areas.  The  York  and  Patuxent  Rivers  were  chosen  to  represent  Virginia  and 
Maryland,  respectively,  because  these  well-studied  river  systems  have  the  necessary  data  to 
calculate  sediment  budgets.  The  study  revealed  problems  unique  to  datasets  for  each  state  during 
the  first  phase  of  the  project.  With  data  issues  resolved,  the  Potomac  River  was  chosen  because  it 
is  the  largest  estuary  of  the  Bay  (the  Bay  is  the  estuary  for  the  Susquehanna  River)  and  the  focus 
of  many  state  and  Eederal  agencies.  Sediment  data  from  the  Chesapeake  Bay  Program  (CBP) 
water  quality  monitoring  stations  in  the  tributaries  were  used  to  study  sediment  transport. 
Stations  were  chosen  by  location  (Eigure  2).  The  “original”  stations  were  the  head  and  mouth  of 
the  estuaries,  and  the  “additional”  stations  were  added  for  a  better  understanding  of  sediment 
transport  within  the  estuaries.  The  water  quality  monitoring  stations  have  CBP  identifications,  in 
which  the  prefix  “EE”  denotes  “lower  estuary”  and  “RET”  denotes  “river-estuary  transition.” 

While  the  watersheds  of  these  rivers  differ  greatly  in  relief,  areal  extent,  soil  type,  land  use,  etc., 
the  estuaries  for  all  the  systems  are  limited  to  the  low-relief  Coastal  Plain.  Eand  use  ranges  from 
predominantly  forested  in  the  York  River  watershed  to  a  mix  of  developed,  agricultural  and 
forested  in  the  Patuxent  River  and  Potomac  River  watersheds.  Water  quality  monitoring  stations 
used  to  calculate  estuarine  loads  were  located  only  in  the  brackish  part  of  each  system.  All  the 
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systems  have  a  main  channel  that  varies  in  depth  and  tends  to  shallow  upstream,  flanked  by 
broad  shoals  about  2  m  deep. 


Figure  1.  Conceptual  sediment  budget  for  Chesapeake  Bay  estuaries. 


Sediment  Data.  In  this  study,  suspended  sediment  was  represented  using  total  suspended 
solids  (TSS),  which  consists  of  both  inorganic  and  organic  particles.  The  US  Geological  Survey 
(USGS)  is  now  using  suspended  sediment  concentrations  (SSC)  at  Fall  Line  monitoring  stations, 
but  these  are  essentially  unavailable  in  the  estuarine  portions  of  the  Bay.  Due  to  laboratory 
processing,  TSS  samples  are  skewed  towards  finer  sediments  compared  to  SSC.  While  they  are 
not  precisely  comparable,  and  there  is  no  easy  conversion  factor  (Gray  et  al.  2000,  Glysson  et  al. 
2000),  TSS  are  measured  Bay-wide  and  many  stations  have  a  long  period  of  record  (tens  of 
years).  Because  of  data  availability  and  for  consistency,  TSS  were  used  everywhere  in  this  study. 

SEDIMENT  BUDGET  COMPONENTS. 

Estuarine  Sediment  Loads.  Calculating  suspended  sediment  loads  in  estuaries  is 
difficult  because  of  the  complex  nature  of  water  movement  in  these  systems.  Multiple  estuarine 
transport  processes,  such  as  gravitational  circulation,  tidal  pumping,  and  river  flow  all  contribute  a 
portion  of  the  suspended  sediment  load.  In  this  text  (for  uniformity),  the  processes  are  called 
gravitational  transport,  tidal  transport,  and  river  transport,  respectively.  The  most  commonly  used 
sediment  loads  for  Virginia  tributaries  (Schubel  and  Carter  1976)  were  based  on  a  salt  balance 
using  a  single  year’s  data  at  a  station  in  the  main  stem  of  the  Bay.  A  new  methodology  using 
sediment  data  from  tributary  water  quality  monitoring  stations  was  developed  subsequently 
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Figure  2.  Estuaries  and  water  quaiity  monitoring  stations.  LE  denotes  iower  estuary; 

RET  denotes  river-estuary  transition. 
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(Herman  2001)  and  improved  on  here  to  estimate  several  important  aspects  of  suspended 
sediments  in  estuaries,  including  the  magnitudes  and  directions  of  estuarine  transport  processes  and 
the  calculation  of  suspended  sediment  loads  (mean  annual  values)  in  the  estuaries. 

The  basic  concept  is  that  concentration  multiplied  by  velocity,  and  then  integrated  over  the  cross- 
sectional  area,  equals  the  sediment  load. 

Q  =  j  (u  *  c  *  dA)  Equation  (1) 


where  Q  is  load,  u  is  velocity,  c  is  concentration,  A  is  cross-sectional  area,  and  dA  represents  the 
components  of  the  total  area  over  which  the  integration  is  summed. 

For  this  study,  it  was  assumed  that  the  tributaries  have  a  typical  two-layer  circulation  in  a 
partially  stratified  estuary  (Haas  1977,  Owen  1969,  Lin  and  Kuo  2001).  In  a  tidal  system  with 
two-layer  circulation,  the  load  for  the  upper  and  lower  layers  must  be  integrated  separately 
before  combining  into  a  total  load.  It  is  assumed  that,  in  each  of  the  two  layers,  it  is  reasonable  to 
represent  the  velocity  and  concentration  via  representative  layer- averaged  values. 


Qlayer  ~  [(Utidal  transport  ^  C)  +  (Ugravitational  transport  ^  C)  +  (Unver  transport  *  c)]  *  A 

layer 


Equation  (2) 


Qtotal  —  Qupper  layer  +  Qlower  layer 


Equation  (3) 


where  u  (m/s)  was  calculated  from  Chesapeake  Bay  hydrodynamic  model  output  (1985-2005),  c 
(mg/1,  converted  to  g/m  )  is  TSS  concentrations  from  Virginia  and  Maryland  water  quality 
monitoring  stations,  and  A  (m^)  was  calculated  using  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  bathymetric  data. 

The  velocities  attributed  to  tidal  and  gravitational  transport  were  extracted  from  the 
hydrodynamic  model  output,  and  river  discharges  from  USGS  gauging  stations  were  used  for  the 
river  transport.  For  each  transport  process,  the  data  from  the  1985-2005  time  period  were  plotted 
on  a  single  graph  versus  tidal  phase.  Figure  3  shows  sample  graphs  for  station  LE4.2  in  the  York 
River.  These  displays  are  useful  for  elucidating  patterns  in  sediment  transport.  For  example, 
multiplying  TSS  concentrations  by  tidal  velocities  yields  the  tidal  transport  graph  that  shows 
greater  transport  on  flood  tide  in  the  lower  layer  of  the  York  River  at  station  LE4.2,  off 
Gloucester  Point. 


Estuarine  Sediment  Accumuiation.  Sediment  accumulation  in  estuaries  was 
calculated  using  bathymetric  data  from  two  distinct  time  periods  that  differ  for  each  river.  For  the 
York  River,  the  dates  were  1857  and  1945;  for  the  Patuxent  River,  the  dates  were  1944  and 
1985;  and  for  the  Potomac  River,  the  dates  were  1862  and  1955.  Bathymetric  data  were  from 
hydrographic  surveys  conducted  by  the  NOAA  National  Ocean  Service  and  are  available  in 
digital  format  online  (National  Oceanic  and  Atmospheric  Administration  2012).  Some  of  the 
older  surveys  are  unavailable  in  digital  format  and  were  digitized  previously  by  the  Virginia 
Institute  of  Marine  Science  staff.  The  soundings  for  the  older  dataset  were  adjusted  for  the 
difference  in  relative  sea  level  rise  with  the  newer  dataset.  The  two  sets  of  surveys  were  each 
converted  into  a  bathymetric  surface  with  ArcGIS,  a  geographic  information  system,  using  an 
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LE4.2  -  TSS  Concentrations  -  Lower  Layer 

n=  131  n=  119 

♦ 

000 

A 

♦ 

♦  ♦ 

u 

♦ 

♦ 

♦ 

♦  .f 

♦ 

^  ♦ 

♦  ♦ 

♦  ♦ 

*  ♦ 

%**♦***>* 

t  ♦  ♦ 

X,-  ’. 

-  V 

s 

0.0 

ack  eb 

DOO  0.2 

30.681  ± 

D  slack  floe 

500  0.5000  0.7 

2.265  Tidal  Phase  57.56 

44.124  ±  3.009 

d  slack 

500  1.0 

7  ±  5.576 

0.4  - 

0.3  - 

0.2  - 

0.1  - 

'w  0  - 

0.0 
3  -0.1  - 

-0.2  - 

-0.3  - 

-0.4  - 

-0.5  - 

LE4.2  -  Tidal  Velocity  -  Lower  Layer 

n=131  n=  119 

»♦ 

)00 

Jr/k 

tvUV 

♦♦ 

7^ 

• 

)00^  ^  0.2 

)00  0.7 

)00  1.0 

w 

‘  y  ♦  ♦ 

♦  ♦♦♦♦ 

slack  ebb  slack  flood  slack 

-0.173  ±  0.008  Tidal  Phase  0.173  ±  0.008 

0.000  ±  0.006 

LE4.2  -  Tidal  Transport  -  Lower  Layer 
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Figure  3.  TSS  concentrations,  tidai  veiocities,  and  tidai  transport  for  station  LE4.2. 
Means  and  standard  errors  of  ebb  and  fiood  phases  are  in  biack.  Totai  mean  and 
standard  errors  are  in  biue.  When  appiying  the  resuits  of  this  study,  it  is 
recommended  that  aii  vaiues  be  rounded  to  two  significant  figures. 
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interpolation  method  that  creates  triangulated  irregular  networks  (TINs).  The  two  surfaces  were 
then  subtracted  and  a  volume  change  (m^)  was  calculated  (Figure  4). 


Figure  4.  Estuarine  sediment  accumuiation  in  the  York  River  using  TINs  created  with 
bathymetry  from  1857  and  1945.  Upper  inset  shows  iocation  of  channeis  (dark  gray) 
and  shoais  (iight  gray).  Lower  inset  shows  iocation  of  bathymetric  change  data. 
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Shoreline  Erosion.  Shoreline  erosion,  SE,  was  calculated  by  determining  the  volume  of 
material  input  along  each  reach  of  undefended  shoreline: 

SE  =  [(reach  length)  *  (average  bank  height)  *  (average  erosion  rate)]  Equation  (4) 

Data  came  from  several  sources  including  historical  erosion  rates  (Maryland  Geological  Society 
2012)  and  EIDAR  (Eight  Detection  And  Ranging)  data  from  Maryland,  and  digital  terrain  models 
(State  of  Virginia  2012)  and  historical  bank  erosion  studies  (Hardaway  et  al.  1992)  from  Virginia. 

Shoreline  erosion  consists  of  fastland  (subaerial)  and  nearshore  (subaqueous)  erosion.  Only 
fastland  erosion  was  calculated  for  the  budgets  because  nearshore  erosion  often  is  a  much 
smaller  proportion  of  shoreline  erosion  (Hardaway  et  al.  2009).  Additionally,  nearshore  erosion 
is  captured  in  part  by  the  estuarine  sediment  accumulation  term,  and  separating  eroding  sediment 
from  resuspended  sediment  in  the  nearshore  is  problematic. 

Biogenic  Production.  Currently,  it  is  thought  that  biogenic  production  may  be  a 
substantial  portion  of  sediment  sources  in  mid-Bay  estuaries  (Cronin  et  al.  2003).  Biogenic  silica 
values  were  collected  from  numerous  samples  in  the  Rappahannock  River  (Anderson  1982). 
Concentrations  for  all  stations  in  the  estuary  for  all  seasons  were  averaged  to  obtain  a  mean 
concentration  of  0.15  mg/1.  The  contribution  from  biogenic  production  using  this  mean 
concentration  is  quite  small;  in  the  Patuxent  River,  it  is  less  than  1%  of  the  smallest  source. 

Another  important  consideration  is  the  use  of  TSS  data,  which  already  incorporates  most  or  all  of 
the  biogenic  production  so  that  a  separate  term  in  the  sediment  budget  was  deemed  unnecessary. 
Both  TSS  and  SSC  contain  inorganic  sediment  (e.g.,  sand,  silt,  clay),  organic  material  from 
runoff  (e.g.,  leaves,  peat)  and  in  situ  biogenic  production  (e.g.,  diatoms,  foraminifera).  Eixed 
suspended  solids  (ESS)  are  the  remains  when  TSS  samples  are  heated  to  remove  the  volatile 
organic  material.  Although  the  volatile  material  is  burned  off,  the  organic  tests  of  the 
microorganisms  remain.  Therefore,  the  biogenic  production  term  in  these  sediment  budgets  was 
assigned  a  value  of  0. 

DISCUSSION.  Sediment  budgets  for  the  York  River,  Patuxent  River,  and  Potomac  River 
(partial  budget)  are  shown  in  Eigures  5  through  7.  To  illustrate  the  dynamic  nature  of  these 
systems,  sediment  budgets  include  intermediary  loads  in  the  estuaries  as  well  as  those  at  the  head 
and  mouth.  The  results  of  this  study  represent  the  most  comprehensive  calculations  to  date  of 
sediment  loads  for  Bay  tributaries. 

Although  the  rivers  show  different  magnitudes  and  transport  directions  of  sediment  loads, 
averaging  all  the  loads  for  each  river  shows  a  net  input  of  sediment  into  each  estuary,  supporting 
the  traditional  thinking  that  tributaries  are  sinks  of  sediment  from  the  Bay  (Eigure  8). 
Additionally,  the  absolute  value  of  the  total  average  load  for  each  river  decreases  from  south  to 
north  in  the  Bay  especially  when  normalized  by  river  cross-sectional  area.  Total  average  loads 
from  the  James  and  Rappahannock  Rivers  are  needed  to  confirm  or  challenge  this  pattern.  This 
tendency  for  decreasing  absolute  load  to  the  north  is  consistent  with  decreasing  tidal  range  along 
this  portion  of  the  Bay  main-stem  and,  thus,  decreasing  magnitude  of  tidal  resuspension. 
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Budgets  for  the  York  and  Patuxent  Rivers  both  show  a  sediment  loss  that  is  unaccounted  for  (i.e., 
to  balance  the  budget  so  the  error  term  is  zero,  more  sediment  is  needed  from  sources,  or  the  sinks 
are  too  large  and  need  to  be  reduced).  Many  factors  contribute  to  the  error  term  in  a  sediment 
budget  mass  balance  equation.  The  main  causes  are  spatial  scales,  temporal  scales,  and  missing 
data.  All  efforts  were  made  to  recognize  these  differences  and  the  effects  they  may  have  on  the 
results.  One  way  to  help  reduce  inconsistencies  is  to  calculate  a  value  for  each  component  of  the 
sediment  budget,  so  that  no  components  were  determined  by  subtraction.  This  puts  the  errors  from 
all  the  components  into  one  term,  which  then  can  be  analyzed  and  apportioned  to  possible  causes. 


Figure  5.  York  River  estuarine  sediment  ioads  and  budget.  Loads  are  in  Mt/yr. 


Long-term  pre-existing  datasets  provide  unique  challenges,  but  also  provide  a  different 
perspective  from  using  short-term  localized  monitoring.  Ideally,  data  for  all  budget  components 
would  be  collected  over  the  same  time  span,  but  in  reality,  this  rarely  occurs.  In  general, 
geomorphic  rates  for  shorter  time  spans  are  higher  and  more  variable,  so  the  longer  time  spans 
may  be  more  comparable.  Despite  these  issues,  the  results  in  this  study  are  based  on  the  best 
available  data.  As  new  data  are  collected,  future  iterations  of  the  budgets  will  improve. 
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Figure  6.  Patuxent  River  estuarine  sediment  ioads  and  budget.  Loads  are  in  Mt/yr. 


MANAGEMENT  APPLICATIONS.  The  sediment  budgets  developed  in  this  study  allow 
an  increased  understanding  of  the  patterns  of  sediment  transport,  erosion,  and  deposition  in 
Chesapeake  Bay  estuaries.  For  example,  in  the  lower  York  River,  the  large  difference  in  loads 
between  stations  LE4.3  and  LE4.2  suggests  a  tendency  for  erosion  between  these  two  stations 
(Eigure  8a).  However,  there  is  not  enough  shoreline  erosion  (most  plausible  long-term  source  of 
sediment)  to  account  for  the  difference.  The  York  River  in  particular  is  highly  energetic,  with 
high  levels  of  resuspended  sediment  (Dellapenna  et  al.  1998)  that  may  account  for  a  substantial 
portion  of  the  calculated  load.  Thus,  bed  sediment  may  be  exported  from  between  these  two 
stations  during  the  non-storm  conditions  captured  by  monitoring. 
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Figure  7.  Potomac  River  estuarine  sediment  ioads  and  partiai  budget.  Loads  are  in 

Mt/yr. 


In  the  Patuxent  River,  the  large  load  at  station  LEI. 2  may  be  due  in  part  to  sediment  input  from 
the  adjacent  tributary.  Sediment  from  upland  erosion  in  the  Coastal  Plain  usually  is  input  through 
streams.  It  is  thought  that  the  significance  of  upland  erosion  in  this  type  of  sediment  budget  is 
minimal  because  of  the  low  relief  (Gellis  et  al.  2007).  The  relative  contribution  of  the  erosion  of 
sediment  from  the  land  surface  rather  than  from  stream  corridors  is  not  well  understood  in  the 
Chesapeake  Bay  basin  (Gellis  et  al.  2007).  Further  investigation  of  sediment  distribution  in  these 
systems  would  contribute  information  about  cumulative  input  from  the  numerous  small 
tributaries  of  an  estuary,  and  help  with  water  quality  and  shoreline  management  decisions. 

The  methods  developed  and  used  in  this  study  are  transferable  to  other  systems.  Sediment  loads 
and  sediment  budgets  from  other  rivers  in  the  Bay  would  help  clarify  the  overall  picture  of 
sediment  transport  and  distribution,  and  assist  regional  sediment  management  efforts  in 
Chesapeake  Bay. 
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Figure  8.  Magnitude  and  direction  of  estuarine  sediment 
loads  for  each  station,  and  of  average  loads  for  each  river 
(flood  is  positive;  ebb  is  negative),  (continued) 
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Figure  8  (continued).  Magnitude  and  direction  of  estuarine 
sediment  loads  for  each  station,  and  of  average  loads  for 
each  river  (flood  is  positive;  ebb  is  negative). 
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ADDITIONAL  INFORMATION.  This  Coastal  and  Hydraulics  Engineering  Technical  Note 
(CHETN)  was  prepared  as  part  of  the  Regional  Sediment  Management  (RSM)  program,  and  was 
written  by  Dr.  Julie  D.  Herman,  Center  for  Coastal  Resources  Management,  Virginia  Institute  of 
Marine  Science  (VIMS),  Gloucester  Point,  VA,  for  the  US  Army  Engineer  District,  Baltimore 
(NAB).  Herman  and  Eriedrichs  (2010)  contains  additional  information,  data,  and  figures,  and  is 
available  at  http://ccrm.vims.edu/research/water  sediments/suspended  sed i me nts/i ndex.htm I 

Carl  Eriedrichs,  Department  of  Physical  Sciences,  VIMS,  collaborated  on  calculating  the 
estuarine  sediment  loads.  This  project  was  supported  by  NAB  as  part  of  the  Corps’  Regional 
Sediment  Management  (RSM)  program.  This  work  supports  the  RSM  efforts  in  the  Chesapeake 
Bay  for  both  the  Baltimore  District  and  Norfolk  District.  The  RSM  point  of  contact  (POC)  for 
the  Baltimore  District  is  Michele  Gomez.  Additional  information  pertaining  to  RSM  can  be 
found  at  the  Regional  Sediment  Management  web  site  http://rsm.usace.armv.mil 

Questions  pertaining  to  this  CHETN  may  be  directed  to: 

Dr.  Julie  D.  Herman  herman@vims.edu 

Michele  Gomez  Michele. Gomez@usace. army. ml 

(RSM  POC  at  NAB) 

Linda  S.  Lillycrop  Linda.S.Lillvcrop@usace.armv.mil 

(RSM  Program  Manager) 

This  ERDC/CHE  CHETN-XIV-23  should  be  cited  as  follows: 

Herman,  J.  D.  2012.  Estuarine  sediment  budgets  for  Chesapeake  Bay  Tributaries. 

Coastal  and  Hydraulics  Engineering  Technical  Note  ERDC/CHE  CHETN-XIV- 
23.  Vicksburg,  MS:  US  Army  Engineer  Research  and  Development  Center, 

Coastal  and  Hydraulics  Eaboratory. 

http://chl.erdc.usace.armv.mil/librarv/publications/chetn/pdf/chetn-xiv-23.pdf 
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ACRONYMS  AND  ABBREVIATIONS. 

Term 

CBP 
CHETN 
CHL 
ERDC 
ESS 
LE 

LIDAR 
NAB 
NOAA 
RET 
RSM 
SSC 
TIN 
TSS 
USAGE 
USGS 
VIMS 


NOTE:  The  contents  of  this  technical  note  are  not  to  be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official  endorsement  or 

approval  of  the  use  of  such  products. 


Definition 

Chesapeake  Bay  Program 

Coastal  and  Hydraulics  Engineering  Technical  Note 

Coastal  and  Hydraulics  Laboratory 

Engineer  Research  and  Development  Center 

Fixed  Suspended  Solids 

Lower  Estuary 

Light  Detection  And  Ranging 

North  Atlantic,  Baltimore 

National  Oceanic  and  Atmospheric  Administration 

River-Estuary  Transition 

Regional  Sediment  Management 

Suspended  Sediment  Concentrations 

Triangulated  Irregular  Networks 

Total  Suspended  Solids 

US  Army  Corps  of  Engineers 

US  Geological  Survey 

Virginia  Institute  of  Marine  Science 
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